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We develop a theory ol the conductance of a quantum dot connected to two leads by single-mode 
quantum point contacts. If the contacts are in the regime of perfect transmission, the conductance 
shows no Coulomb blockade oscillations as a function of the gate voltage. In the presence of small 
reflection in both contacts, the conductance develops small Coulomb blockade oscillations. As the 
temperature of the system is lowered, the amplitude of the oscillations grows, and eventually sharp 
periodic peaks in conductance are formed. Away from the centers of the peaks the conductance 
vanishes at low temperatures as T 2 , in agreement with the theory of inelastic co-tunneling developed 
for the weak-tunneling case. Conductance near the center of a peak can be studied using an analogy 
with the multichannel Kondo problem. In the case of symmetric barriers, the peak conductance at 
T — > is of the order of e 2 /%. In the asymmetric case, the peak conductance vanishes linearly in 
temperature. 



PACS numbers: 73.20.Dx, 73.40.Gk 

I. INTRODUCTION 

Electronic transport in mesoscopic systems is usually 
affected by the Coulomb interactions between the elec- 
trons. The interactions manifest themselves most dra- 
matically in small tunnel junctions. If a small metallic 
grain is connected to a lead by a tunnel junction, an 
electron tunneling into such a grain charges it by the 
elementary charge e and increases the energy of the sys- 
tem by e 2 /2C. For a small grain, the typical capaci- 
tance of the system can be very small, so that the charg- 
ing energy e 2 /2C can significantly exceed the tempera- 
ture of the system T. In this regime one observes the 
Coulomb blockade: strong suppression of the tunneling 
into the grain at low temperatures. Various aspects of 
the Coulomb blockade have been intensively studied, in 
recent years, both experimentally ancLtheoreticallyoEl 

Recently it has become possible^ to observe the 
Coulomb blockade in semiconductor heterostructures. 
Instead of a metallic grain, in such experiments one cre- 
ates a small region of the two-dimensional electron gas 
(2DEG) — a quantum dot — separated from the rest of the 
2DEG by potential barriers created by applying negative 
voltage to the gates, Fig. |l|. Experimentally, one mea- 
sures the linear conductance through the quantum dot as 
a function of the voltage V g applied to the central gates. 
The role of the gate voltage is to change the electrostatic 
energy of the system, 

(Q_eNl 

where Q is the charge of the quantum dot, and N is a di- 
mensionless parameter proportional to the gate voltage 
and the capacitance C g between the dot and the gate, 
N = CgVg/e. If the potential barriers separating the 
quantum dot from the leads are high, the charge Q is 
integer in units of e. Thus, the ground state of the sys- 
tem corresponds to some integer value of charge, and the 



states with different values of Q are separated from it 
by the electrostatic gap A ~ e 2 /2C. It is important to 
note that at the values of the gate voltage correspond- 
ing to half- integer values of N the gap vanishes, i.e., the 
ground state is degenerate in charge Q. For example, at 
N = \ the ground state can have either charge Q = or 
e. Therefore at half- integer N the tunneling of an elec- 
tron into or out of the quantum dot does not lead to the 
increase of the electrostatic energy of the system, and 
the Coulomb blockade is lifted. As a result the conduc- 
tance of the structure in Fig. |l| shows periodic peaks as 
a function of N. The Coulomb blockade peaks in linear 
conductance as a function of the gate voltage are readily 
observed in the experiments (see, e.g., Ref. @). 




FIG. 1. Schematic view of a quantum dot connected to two 
bulk 2D electrodes. The dot is formed by applying negative 
voltage to the gates (shaded). Solid line shows the boundary 
of the 2D electron gas (2DEG). Electrostatic conditions in the 
dot are controlled by the voltage applied to the central gates. 
Voltage Vl,r applied to the auxiliary gates controls the trans- 
mission probability through the left and right constrictions. 
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The shapes and heights of the peaks in conductance 
are usually described within the framework of the ortho- 
dox theory of the Coulomb blockaded In this approach 
the conductance is found from the balance of the rates 
of tunneling of electrons between the leads and the dot, 
with the rates being calculated in the lowest order per- 
turbation theory in the tunneling matrix dements. The 
resulting linear conductance has the form:Q 



q _ 1 GlGr 



A/T 



2G L +G R sinh(A/T)' 



A = -5N. 



(2) 



Here Gl and Gr are the conductances of the left and 
right barriers, 6N is the deviation of N from the nearest 
half-integer value. Away from the center of the peak, 
the tunneling of an electron into the quantum dot is 
suppressed, because only a small fraction of electrons 
has the energy sufficient to overcome the electrostatic 
gap A. Thus the conductance (||) decays exponentially 
at T — ► 0. At very low temperatures, however, an- 
other mechanism — usually referred to as inelastic co- 
tunneling — dominates the conductance. Instead of a real 
process of tunneling from a lead to the dot, electron tun- 
nels from the left lead to a virtual state in the dot, and 
then another electron tunnels from the dot to the right 
lead. As a result of such a two-step process, one electron 
is transferred from the left lead to the right one, and 
the charge of the quantum dot is unchanged. The corre- 
sponding conductance has only a power-law dependence 
on temperature,u 



G = 



ttTiGlG! 
3e 2 



T< A. 



(3) 



In the case of weak tunneling, Gl.Gr <§; e 2 /fr, the re- 
sults (|J) and (||) give the full description of the Coulomb 
blockade peaks in linear conductance. Near the center of 
the peak the conductance is given by Eq. (||), whereas 
the tails of the peaks are dominated by the co-tunneling 
contribution (||). n n 

In a number of recent experimentsHij the tunneling 
through quantum dots was studied in the regime when 
one or more quantum point contacts were tuned to the 
regime of strong tunneling, corresponding to conduc- 
tance G — > e 2 /irh. Inj-particular, it was demonstrated 
by van der Vaart et alB that the Coulomb blockade os- 
cillations of linear conductance persist as long as the 
conductances of both barriers are below the conduc- 
tance quantum e 2 /irh. Most of the existing theoretical 
works OB-.the strong-tunneling regime of the Coulomb 
blockadaldilj discuss the case of metallic systems, where 
one can achieve large conductance of a tunnel junction 
between the grain and a lead by increasing the number of 
transverse modes while keeping the transmission coeffi- 
cient T small. Such theories cannot be applied to the ex- 
periments in semiconductor heterostructures, where the 
point contacts typically allow only a single transverse 
mode, but the transmission coefficient can be tuned to 



T — > 1. The theoretical works for this case either da 
not concentrate on the Coulomb blockade oscillationstfl 
or studv-jthe oscillations of equilibrium properties of the 
system,E_3 rather than its conductance. 

In this paper we develop a theory of the Coulomb 
blockade oscillations of conductance in the regime of 
strong tunneling. It was shown in Rcf. ^| that the av- 
erage charge of a quantum dot connected to a lead by a 
single-mode quantum point contact shows small periodic 
oscillations as a function of the gate voltage as long as 
the transmission coefficient T is below unity. Similarly, 
we will show that the conductance through a quantum 
dot shows small periodic oscillations at T < 1. Opa 
should note, however, that unlike the average charge,li3 
the amplitude of the oscillations of conductance depends 
on temperature and grows as T is lowered. In the limit 
T — ► 0, small oscillations develop into sharp periodic 
peaks, and the off-peak conductance vanishes as T 2 , thus 
reproducing the characteristic quadratic temperature de- 
pendence for inelastic co-tunneling 

In section [n] we use the analogyE3EZl between the 
Coulomb blockade and multichannel Kondo problem to 
discuss the peak value of conductance at T — > 0. This 
discussion is done for the weak-tunneling case, but the 
conclusion that the peak conductance is of the order of 
e 2 /h is valid for the strong tunneling as well. The the- 
ory of conductance oscillations in the strong-tunneling 
regime is developed in sections III and [V . In section |v| 
we discuss the scaling properties of our problem. In par- 
ticular, we show that the quadratic temperature depen- 
dence of the off-peak conductance at T — > is universal 
and holds for arbitrary transmission coefficients. 



II. CONDUCTANCE PEAKS IN THE 
WEAK-TUNNELING CASE 



Let us first consider the weak-tunneling case. The lin- 
ear conductance is given by Eqs. (0) and One should 
note that in the derivation of Eqs. (||) and (0) only the 
lowest-order terms of the perturbation theory in tunnel- 
ing amplitudes were taken into account. This is usually 
justified at Gl,Gr <C e 2 /h. However, it was shown in 
Ref. [if] that near a half-integer TV and at low tempera- 
tures the higher-order corrections lead to large logarith- 
mic renormalizations of the tunneling amplitudes. The 
nature of these renoijiaalizations can be understood in 
terms of the analogyO between the Coulomb blockade 
and the Kondo problem. We will now discuss this anal- 
ogy and correct the results (||) and (|J) to take into ac- 
count the renormalizations of the tunneling amplitudes. 

For simplicity we start with the quantum dot weakly 
coupled to a single lead. Such a system can be described 
by the following Hamiltonian: 



H 
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(4) 



Hero a k and a p are the annihilation operators for the elec- 
trons in the lead and in the dot respectively, and Ec is 
the characteristic charging energy, Ec — e 2 /2C; coupling 
between the dot and the lead is described by tunneling 
Hamiltonian with the matrix element v t . Finally, n is the 
operator of the number of electrons in the dot, counted 
from that at the ground state of the system without tun- 
neling, n = Y,( a t a P ~ (4s)o)- 

Due to the last term in Eq. ([!]) electrons can tunnel 
from the lead to a virtual state in the dot. The energy of 
such a virtual state is of the order of Ec, and to complete 
this scattering process an electron has to tunnel from the 
dot to the lead. It is instructive to calculate the second- 
order amplitude of scattering between two states in the 
lead, k and k' , near the Fermi level: 



v ■■ 
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\v t \ 2 9(- 
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(5) 



Here the second term corresponds to the process in which 
first an electron from the filled state p in the dot tunnels 
to the state k' in the lead, and then the electron k fills 
the hole in the dot. Clearly, the total amplitude diverges 
logarithmically at N 



T, 



k^k' 



- v d\v t 



! ln- 



1 + 27V 
1 - 27V' 



(6) 



where v& is the density of states in the dot. The logarith- 
mic singularity at N — > g originates from the degeneracy 
of the electrostatic energy for states with n = and n = 1 
extra electrons in the dot. 

Let us now concentrate on the case of the dimension- 
less gate voltage N very close to \. In this case the 
logarithmic renormalizations are large, and to find the 
scaling properties of the scattering amplitudes, we can re- 
strict our consideration to the energy scales smaller than 
Ec- Thus when constructing the perturbation series for 
a tunneling amplitude we will neglect all the processes 
in which the system visits virtual states with the charge 
n different from or 1. Another important property of 
the Hamiltonian (Q) is that the perturbation — the tun- 
neling term — always changes the charge of the dot by 1. 
Suppose that at the first step of a high-order tunneling 
process an electron tunnels from the lead to the dot thus 
changing n from to 1. Then at the next step an electron 
will have to tunnel from the dot to the lead (changing 
n = 1 — ► 0), because otherwise a virtual state with n = 2 
and correspondingly large energy ~ Ec would have been 
created. At the third step an electron tunnels from the 
lead to the dot, changing n — — > 1, and so on. 

One can easily see that exactly the same perturba- 
tion series for scattering amplitudes is obtained from the 
Hamiltonian of the anisotropic Kondo model, 



H = tka\ a a ka + hS z 



a ka a k>a>i 



(7) 



kk' 



where a ± are the linear combinations of Pauli matrices, 
a = <T x ±a v , and S is the impurity spin, S — S X ±S V . 
A formal discussion of this mapping can be found in 
Rcf. [I?]; here we only mention the relations between the 
parameters of the Hamiltonians (Q) and (0). The values 
of the electron spin a = | and | correspond to the elec- 
tron being in the lead and the dot respectively. Similarly, 
the values of the impurity spin S z = J. and | correspond 
to the states with the number of particles in the dot n = 
and 1 . Magnetic field h is identified with the energy split- 
ting of the n = and n = 1 states, h = 2Ec{\ ~ N). 
Finally, the tunneling matrix element v t is mapped to the 
exchange constant Jj_. 

It is well known that at low temperatures and mag- 
netic fields the exchange constant J of the Kondo model 
renormalizes logarithmically. As a result J grows, and 
the system approaches the strong-coupling fixed point. 
In the context of the Coulomb blockade problem, this 
means that at N — > | and T — > the tunneling ampli- 
tudes grow until the conductance of the barrier becomes 
of order e 2 /h. 

Using the analogy with the Kondo problem, one can 
actually find the renormalizations of the tunneling ampli- 
tudes. In the leading logarithm approximation oue, can 
use the poor man's scaling technique£3 The resulttL£l can 
be written down in terms of the renormalization of the 
transmission coefficient T of the barrier between the lead 
and the dot, 



T = 



=2/1 



e = ln 



Er 



max{T, A} ' 



(8) 



Here T is the bare value of the transmission coefficient, 
and A = Ec\l — 2N\. When T becomes of the order of 
unity, the leading logarithm approximation fails. Thus 
the result (g) is valid as long as the renormalized trans- 
mission coefficient is small, T <C 1. 

In this regime one can still use the weak-tunneling re- 
sults (g) and (||) for the conductance through a dot con- 
nected to two leads, taking into account the renormal- 
izations of the conductances of the barriers. This can be 
done by noting the relation between the conductance of 
a barrier Gb and its transmission coefficient, 



G b = — x 
2nn 



r, 

2T, 



without spins, 
with spins. 



(9) 



It is worth mentioning that within the leading logarithm 
approximation the transmission coefficients of two barri- 
ers renormalize independently, so one can use Eq. (|^) for 
both left and right barriers. 
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Let us find the temperature dependence of the peak 
value of the conductance in the case of symmetric barri- 
ers, i.e., assume N = | and Gl — Gr = Gb — (e 2 /n?i)T. 
According to Eq. (^), the peak conductance is |G&. Tak- 
ing into account the renormalization (||) of the transmis- 
sion coefficient, we find 



G = 



To 



4nh cos 2 [±y/%ln(Ec/T)] 



(10) 



As expected, when the temperature is lowered, the con- 
ductance grows. This result is valid only at temperatures 
exceeding the characteristic Kondo temperature of this 
problem, Tk — Ec exp(— 7r 2 / 2 \/7q). At too low temper- 
atures, T < Tk, the conductance is of the order of e 2 /Ti, 
transmission coefficient T ~ 1 , and the result of the lead- 
ing logarithm approximation (jl^) is no longer applicable. 

It is also interesting to study the tails of the conduc- 
tance peaks. If the gate voltage N is not precisely half- 
integer, then at low temperatures, T«A = 2Ec\N— 1|, 
one can still use the inelastic co-tunneling result (|J) with 
the values of conductances of the barriers renormalized 
according to Eqs. (||) and (j|). The resulting conductance 
has the form: 



G = 



e 2 T L T R 



T 2 



Zirh A 2 cos 2 [±V7t> cos 2 [i- 



(11) 



Similarly to Eq. (10), the result ( pi] ) is valid only as long 
as T <C 1. In this case the condition allows only the 
values of the gate voltage not too close to the resonance, 
A » T K . 

As we have seen, if one starts with the weak-tunneling 
case, in the vicinity of the resonances in G(V g ) the tun- 
neling amplitude grows, and the system approaches the 
strong-coupling regime. Therefore one can use the re- 
sults of the weak-tunneling approach to find the univer- 
sal features of the peaks in conductance which should 
also persist in the strong-tunneling case. Such features 
are: (i) the T — > value of the peak conductance JlC| ) 
for symmetric barriers is of the order of e 2 /Ti, and (ii) 
away from the centers of peaks the temperature depen- 
dence of the conductance at T — > is G cx T 2 . It is 
worth noting that unlike the temperature dependence of 
the co-tunneling conductance (|^), its 1/A 2 dependence 
of the gate voltage is not universal, see Eq. (p"l|). 



III. STRONG CO-TUNNELING OF SPINLESS 
ELECTRONS 

We will start our treatment of the strong-tunneling 
case with a simplified model of spinless electrons. As 
we will see, in this case a complete solution of the prob- 
lem is possible. Experimentally, the spinless case can 
be achieved by applying an in-plane magnetic field, suffi- 
ciently strong to polarize the electron gas in the constric- 
tions. 



A. Theoretical model 



We assume that each of the quantum point contacts 
in Fig. |l| allows the transmission of electrons in only one 
mode. The transmission coefficients of the two contacts 
Tj, and Tr are assumed to be close to unity. This en- 
ables us to start with the Hamiltonian corresponding to 
perfect transmission, = T R = 1, and then construct 
the perturbation theory in small reflection amplitudes, 
and tr. Following Ref. |lj, we will describe each of the 
quantum point contacts by a system of one-dimenskmal 
(ID) electrons. The formal derivation of this modella is 
based on the fact that if the boundaries of the channel are 
very smooth, the transport thoughJt can be described in 
terms of adiabatic wavefunctionslij having essentially a 
ID form. 

Another simplification of the problem is due to the fact 
that the typical energies involved in the Coulomb block- 
ade phenomenon are of the order of the charging energy 
Ec, which is much smaller than the Fermi energy. This 
allows us to linearize the spectra of the ID electrons near 
the Fermi points and to formulate the model in terms of 
the left- and right-moving fermions: 



Ho = ihvF 



Wh\ (s) V^li (x) - V4, 2 (x)VipL2 (x) 



+ip R1 {x)VipRi{x)-iP R2 {x)VijjR2{x) dx. (12) 



Here and ipL2 are the annihilation operators for the 
electrons of the left contact moving to the left and right 
respectively; ipm and ipR2 describe the left- and right- 
moving electrons of the right contact; V = -§^. We asso- 
ciate the centers of the constrictions with points x = 0, 
meaning that electrons in the dot are described by 4'li(2) 
at x > and by ipm(2) at x < 0. Thus the charging en- 
ergy can be presented as 



H c = E c (h-N) 2 



where 



h = 



- 3 = 1,2 



(13) 



dx. 



Here : X: denotes the normal-ordered operator X. Fi- 
nally, the presence of weak backscattering in the contacts, 
described by small reflection amplitudes |ri|,|rfl| -C 1, 
gives rise to the term mixing the left- and right-moving 
particles, 



H' = hv f 



M4i(0)lM0) + MV4u(0)^ 2 (0) + h.c. 



(14) 
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The Hamiltonian H = Hq+Hc+H' gives the complete 
description of the transport through a quantum dot con- 
nected to two leads. In our model the two contacts are de- 
scribed by two independent ID systems. Thus we neglect 
the possibility of coherent transport of electrons from one 
quantum point contact to the other. Such processes lead 
to the additional mechanism of transport .through the dot 
called elastic co-tunneling. It is known,B however, that 
at temperatures exceeding the level spacing e in the dot 
the elastic co-tunneling contribution to the total conduc- 
tance is much smaller than that of inelastic co-tunneling. 
Thus we will concentrate on the case T ^S> £ which is easy 
to satisfy in the experiments, and in the following calcu- 
lations we will assume e = 0. At T < e the conductance 
becomes temperature-independent and can be estimated 
by replacing T by e in our results. 

The solution of the problem with the Hamiltonian 
given by Eqs. (|l2])-(|l4|) is not trivial because the inter- 
action term Hq is not quadratic in fermion operators. 
To overcome this difficulty, we will use the bosonization 
techniqueEJ which enables one to treat the interaction in 
ID systems exactly. Following Ref. ^0], we present the 
fermion operators as 



Mx) = (^) 1/ % L1 e^ 1 ^, 

1/2 



1pL2{x) 



D 

2-kTivf 



r) L2 eif>L2(x), 



(15) 



Here D is the bandwidth, and the two independent 
bosonic fields tp Ll and (p L2 satisfy the following commu- 
tation relations: 

[<Pli( x )> <PLi(y)] = -[<fL2( x ), ( PL2(y)} = -insgn(x - y). 

(16) 

Given (filf), the operators ( pHl) satisfy usual anticommu- 
tation relations. To ensure proper anticommutation re- 
lations for fermions at different branches, we have added 
the local Majorana fermions rj L1 and rj L2 . 

The electrons of the right contact are bosonized in a 
similar way. One can then rewrite the Hamiltonian (|l^ ) 
( |l4| ) in terms of (fL and ifR. It will be more convenient, 
however, to present the Hamiltonian in terms of slightly 
different bosonic fields: 

<t>L = ~(<PL2- <Pli), TL l = —^-V(<pli + <Pl2)- (17) 

The new variables satisfy the commutation relations typ- 
ical for the displacement of a ID elastic medium and its 
momentum density, 



[c/ )L (x),U L (y)}=td(x-y). 



(18) 



The form of the free-electron Hamiltonian (12) supports 
this observation, 



This qualitative interpretation allows one to predict 
the correct bosonized form of the charging energy jl5|). 
One should expect that the charge of the dot is pro- 
portional to the difference of the charge Ql ~ e<^z,(0) 
brought into the dot through the left contact and the 
charge Qr ~ e<j)R (0) carried away through the right con- 
tact. Explicit calculation gives 



#c = [<M0) - <M0) - TrAf 



(20) 



It is important to stress that the interaction term ( |20| ) 
is now quadratic in 0^ and <j)R. Therefore the bosoniza- 
tion approach enables one to treat the charging energy 
exactly. 

Finally, the backscattering term (R4h takes the form 



H' = cos[2fo(0)] + cos[2fo(0)]. (21) 



In Eq. ( |21| ) we neglected the products of Majorana 
fermions T) L1 rj L2 and tIriVrzi as they commute with each 
other. Unlike the first two terms, H and He, the 
backscattering term is not quadratic in bosonic variables 
and cannot be treated exactly. Below we assume weak 
backscattering and develop the perturbation theory in 
small parameters and |t\r|. 



B. Perturbation theory for the conductance 

In this paper we are only interested in the linear con- 
ductance through the quantum dot. To find it, we will 
use a Kubo formula, which in the zero-frequency limit 
can be written as 



G= 2hl W)> J (0)]>*, 



(22) 



where I is the operator of current through the dot. 

We define I as the average of the current II flow- 
ing into the dot through the left contact and the cur- 
rent Ir flowing out of the dot through the right contact. 
Clearly, the current 7j, is proportional to the time deriva- 
tive of the displacement of the effective clastic medium 
describing the left contact, II ~ e< pL, and, similarly, 
Ir, ~ e(j)R. From the equations of motion correspond- 
ing to the Hamiltonian Hq + He + H', one concludes 
that I ~ evpiJ^L + Hr). A careful calculation gives 



1 = \ ev F (n L + Uj 



lx=0 ' 



(23) 



hvp 



i 



y0i(a:)] 2 +7rll£(a:) 



+ -[V4>r(x)Y + nIl R (x) y dx. 



(19) 



As the first step, we will find the conductance in 
the absence of the barriers. Since the current operator 
( p3| ) commutes with He, one should expect the conduc- 
tance to be insensitive to the charging energy. Clearly, 



5 



at Ec — the system is equivalent to two resistances 
2Trh/e 2 of the two quantum point contacts, connected in 
series. One then expects the conductance to be 



(24) 



This result is easily derived from the Kubo formula (22), 
because with H ' = the current-current correlator must 
be found for a quadratic Hamiltonian. As a result, one 
finds ([!(*), J[0)]) = i(e 2 /2Tr)S'(tV and after the substi- 
tution into ( B2|) o ne reproduces (E4). 

The result (p4|) is independent of temperature. Since 
we used the Hamiltonian (|l^) with linearized spectra 
of electrons, the temperature was assumed to be much 
smaller than the Fermi energy, but could still be of the or- 
der of the charging energy. In the presence of the barriers 
the conductance must acquire some temperature depen- 
dence, which is obvious from the limit of very high bar- 
riers (||). Indeed, the lowest non- vanishing (2nd-order) 
correction in and r/? to the conductance is already 
temperature dependent. At low temperatures T Ec, 
we get 



G = 



ttT (N) 



4T 



(25) 



where T is defined as follows: 



r (N) 



2 1 E C 



[K| 2 + M 2 + 2|r L |Mcos(27rA0]. 



(26) 



Here 7 = e , with C w 0.5772 . . . being the Eulcr's con- 
stant. In Eq. ( p5| ) we have neglected the terms small in 
T /Ec- The result (25) and ( Pq ) is analogous to-that for 
the conductance of a ID wire with two defecta3E3 and 



can be obtained in a similar way, see Appendix A 1 



As expected, the presence of weak scattering in the 
contacts gives rise to a small periodic correction to the 
conductance. The same behavior was predicted^ for the 
average charge of the dot. However, unlike the average 
charge, the correction to conductance strongly depends 
on temperature. As the temperature is lowered, the con- 
ductance decreases. The only case when conductance is 
temperature-independent is when the barriers are sym- 
metric, |r/| = \tr\, and the dimensionless gate voltage 
N is half-integer, which corresponds to the center of a 
peak in the high-barrier case. In this regime we expect 
the conductance to be of the order of e 2 /h, cf. Sec. ||, 
which is confirmed by the perturbation theory ( |25| ) and 



C. Non-perturbative calculation 

At low temperatures T < Tq the second-order correc- 
tion becomes large, and the perturbation theory fails. To 



find the conductance in this case, one has to sum up in- 
finite number of terms of the perturbation theory. To 
perform such a summation, we will use the technique de- 
veloped in Ref. [la. 

At the first step we perform a linear transformation of 
the bosonic fields, 



'L = 



V2 



h + <t>c 
V2 ' 



(27) 



and similarly for and Hr. The advantage of using 
these new variables is that the Coulomb energy fl2C| ) is 
now expressed in terms of the charging mode <pc only, 



Hr 



2E C 



irN 

71 



(28) 



The charging energy ( |28| ) suppresses the fluctuations of 
</>c(0) at energy scales below Ec- Thus at T <C Ec one 
can integrate out the fluctuations of the charging mode 
by replacing the scattering term ( pi] ) with its value av- 
eraged over the fluctuations of (f>c- The resulting Hamil- 
tonian is expressed in terms of the field <fii only and has 
the form 



H' = 



-fE c D 



V2 r 



s jV20/(O) , r * e -n/2<£j(0) 



(29) 



(30) 



where the dependence on the gate voltage is incorporated 
into the complex parameter r = \rL\e~ l7rN + \rn\e l7vN . 

The Hamiltonian (29) and (ROT) can be transformed to 
that of an impurity in a g = ^ Luttinger liquidEa by a 

simple transformation <fi = fJ^4>i and n = II// \/2. The 
latter can be solved exactly^ using the technique devel- 
oped in Ref. ^4[ For our purposes it will be-,more conve- 
nient to use an alternative exact solutionis The idea is 
to interpret e 1 ^ 2 ^ 1 as^bosonized form of some fermion 
annihilation operators This can be done by rewriting 
the Hamiltonian in terms of two decoupled chiral bosonic 
fields defined as 



(p±(x) 



(j>i{x) ± 4>i{-x) 
V2 



n I (x')±ii I (-x') 

V2 



dx'. 



(31) 



In these variables the Hamiltonian takes the form 

H ° = l£j {[ V ^)] 2 + [ V M^)] 2 } dx, (32) 



H' = 



{ 2tt 3 



1/2 



i V+ (0) 



+ r e 



(33) 



One can easily check that the fields <p+ and </?_ commute. 
Thus the </?_ part of the Hamiltonian completely decou- 
ples. Since the field (p+(x) satisfies the same commuta- 
tion relations ( |l6[ ) as <pli, one can in complete analogy 
with Eq. (|l5|) interpret e lip+ as a fermion operator,c3 



G 
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FIG. 2. Conductance (^8|) as a function of the dimension- 
less gate voltage N for the symmetric case, \tl\ = \rn\ = 0.2. 
The two curves are calculated for Ec /T = 25 (a) and 200 (b). 
In this calculat ion we included the correction proportional to 
T/Ec in Eq. (A9). The dashed line represents the conduc- 
tance G = GlGr/(Gl + Gr) in the high-temperature limit, 
T > E c . 



(^fV^). 04) 



The resulting Hamiltonian is quadratic in fermion op- 
erators, 



H = ihvF / ip + (x)V'ip + (x)dx 



n? + ^ + (0) + rVt(0)»7, 



(35) 



and can be easily diagonalizedS This greatly simplifies 
the calculation of the conductance. One can again use 
the Kubo formul a (22 ), with the properly transformed 
current operator (|23 ) . In the new fermionic variables the 
current operator has the form 



I = —cvf 



x=0 



(36) 



The actual calculation of the linear conductance is 
rather long, but straightforward. First we find the 
current-current correlation function to be 



[I(tll(0)])=i-S'(t) 



T 2 sm[(e + e')t/h]dede' 



(37) 



2tt 2 ?L 2 J J (e 2 + r§)(e«/ T + \){e^ T + 1) ' 

The substitution of Eq. ( J37| ) into the Kubo formula (|2^ 
gives the following expression for the linear conductance 



G = 



Airk 



1 - 



dE- 



1 



rg(J\0 



4Tcosh 



2 E 
2T 



E 2 +T 2 {N) 



(38) 



Equation (|38|) is the central result of this section. It 
is valid for any temperature below Ec, and, unlike the 
perturbative result (p5f), it is not restricted to only suffi- 
ciently high temperatures T Tq. Of course, at T 3> To 
the result ( p5| ) is reproduced. 

In the opposite limit, T C To, the conductance is 
quadratic in temperature, 



G 



7re 
12ft 



T 



T {N) 



(39) 



It is important to note that the same temperature de- 
pendence was found for the inelastic co-tunneling (j^) in 
the weak-tunneling regime. As expected, the quadratic 
temperature dependence is a universal property of inelas- 
tic co-tunneling and is not limited to the weak-tunneling 
limit, cf. Eq. (0). 

In the special case of symmetric barriers, \r^ \ = \r^\ = 
ro, the low- temperature limit ( |39| ) is never achieved on 
resonance. Indeed, we have r = (8^ Eq /TT 2 )r 2 cos 2 nN , 
i.e., Tq vanishes at half- integer values of N. Therefore 
in the low-temperature limit the resonant value of con- 
ductance is e 2 /4irh, whereas away from the resonance 
To 7^ 0, and G oc T 2 . It is clear then that at low temper- 
atures the weak Coulomb blockade oscillation in conduc- 
tance ( p5| ) transforms into sharp periodic peaks, Fig. ||. 
Unlike the weak-tunneling case ||), the tails of the peak 
at iV = i are not Lorentzian: G oc 1/(N — i) 4 . 

Finally, it is worth noting that the conductance d3§| ) 
coincides withJhat of a weak impurity in the g = h Lut- 
tinger liquid.ES This is not surprising, considering the 
aforementioned identity of the Hamiltonian (53), ( |30| ) 
with the one for the quantum impurity problemH In the 
latter problem the condition g = \ requires special fine- 
tuning of the strength of interaction between electrons. 
In our case the condition g — \ is satisfied automatically 
when we integrate out one of the two bosonic modes.Ell 



IV. CO-TUNNELING OF ELECTRONS WITH 
SPINS 



In the previous section we considered the case of spin- 
less electrons, which can be realized by applying strong 
in-plane magnetic field. In the absence of the field, each 
of the contacts has two identical modes, corresponding 
to two possible values of electron spin. The number of 
modes is an important parameter in our problem. This 
can be easily seen from the analogy to the Kondo prob- 
lem discussed in Sec. |l[ In the weak-tunneling case, the 
model of spinless electrons maps to the 2-channel Kondo 
problem, because the fictitious spin of the impurity — i.e., 
the charge of the dot — can be changed by tunneling of 
an electron in two independent channels, corresponding 
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to the left and right contacts. In the presence of real 
spins of electrons, the number of channels in the effective 
Kondo problem is 4. It is well known that the num- 
ber of chanacls in the Kondo problem strongly affects its 
properties E3 Thus one should expect that the presence 
of spins should have a dramatic effect upon the transport 
through a quantum dot. 



A. Strong-tunneling case 

We will start with the discussion of the strong- 
tunneling limit, Tl,Tr — > 1. To find the conductance, 



one can use the same technique as in Sec. III. The Hamil 



tonian of the problem can be written in the bosonized 
form, analogous to (pj|)-(pl|), 



;[V0 LQ (a:)] 2 +^nL(x) 



+ -[V<fe a (x)] 2 + Trills) \dx, (40) 

7T 



He 




<j> La (0)]-irN 



(41) 



H ' = ^£{kL|cos[20 ia (O)] + |r fl Jcos[2^, Q (O)]}. (42) 



Here the fields with a = | and J. describe the electrons 
with spins up and down respectively. 

One can again find the linear conductance using the 
Kubo formula (^2j), but the current operator (|23|) must 
be modified to account for the spins, 



(n La + n Ra ) 



(43) 



In the absence of barriers, \tl\ = \tr\ = 0, the cal- 
culation is trivial, as the Hamiltonian ([hj]) and (|4l]) is 
quadratic. Due to the doubling of the number of modes 
in each contact, in the presence of spins the conductance 
doubles compared to Eq. (p4[), 



Go 



e 

2^h' 



(44) 



Weak backscattering, \ri\, \tr\ <C 1, ca n be treated in 
perturbation theory in H' , see Appendix |A 2] , Up to the 
second order we get 



G = 



2-kTi 



1 - 



2£(|) /jEc 



ttT 



(\r L \ 



(45) 



It is instructive to compare this result with its analog 
for the spinless case, namely Eq. (|25|). In both cases the 
small correction to the conductance has the tendency to 
grow as the temperature is decreased. An important dif- 
ference is that the second-order correction in (f45|) does 



not depend on the gate voltage N, and as a result it 
grows even at half-integer N, i.e., on resonance. Thus, 
in contrast to the spinless case, one should expect the 
resonance value to be smaller than the unperturbed con- 
ductance Go = e 2 /2n~h. 

This difference has a natural interpretation in terms 
of the analogy with the Kondo model. As we already 
mentioned, the spinless case corresponds to the 2-channcl 
Kondo model. In the 2-channel case the strong-coupling 
fixed point is stable at zero magnetic field, or, in terms of 
the Coulomb blockade problem, at half-integer N. There- 
fore, precisely on resonance a weak reflection must be an 
irrelevant perturbation, and the peak value of the con- 
ductance should be given by Eq. (|i|). On the contrary, 
the 4-channel Kondo problem corresponding to the case 
with spins has a stable fixed point at an intermediate 
value of coupling, and both weak- and strong-coupling 
fixed points are unstable.cll Thus even on resonance the 
conductance is not given by its value in the strong- 
tunneling limit (44), as confirmed by the perturbation 
theory (Bs). 



B. Conductance in the case of asymmetric barriers 

Unfortunately, it is not easy to calculate the conduc- 
tance near the intermediate-coupling fixed point. Some 
conclusions about it can be drawn on the basis of the 
scaling approach, Sec. [v|. One should note, however, 
that this fixed point is stable only if the coupling in all 
the 4 channels is identical. In the case of the Coulomb 
blockade problem, the coupling in the pairs of channels 
corresponding to different spin orientations in the same 
contact must be identical (in the absence of magnetic 
field) . On the other hand, the coupling constants for the 
modes in different contacts are determined by the corre- 
sponding conductances and are not necessarily equal. In 
fact, in the structure shown in Fig. |l|, the two conduc- 
tances are controlled separately by the gate voltages Vl 
and Vr. Therefore in experiments the coupling constants 
corresponding to the left and right modes are not equal, 
unless a special attempt to make them equal is made. 

In this section we study the limit of very asymmetric 
barriers, namely we assume that one of the barriers is 
very high, whereas the other one is low. In terms of the 
transmission coefficients of the barriers this limit corre- 
sponds to Tl -C 1 and 1 — Tr -C 1. Such a regime can 
be easily realized in an experiment by proper tuning of 
the gate voltages Vl and Vr. In addition, we will see in 
Sec. [V| that at T — ► arbitrarily weak asymmetry grows 
and the system scales towards the strongly asymmetric 
limit. Therefore at low temperatures the results of this 
section can be applied to any asymmetric system. 

In Sec. we described the contact in the regime of 
weak tunneling by the tunnel Hamiltonian (Q) written 
in ter ms o f the original electrons. On the other hand, 
in Sec. Ill it was more convenient to use the bosonized 
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form of the Hamiltonian for the electron gas in the con- 
tact with transmission coefficient close to unity. There- 
fore when considering the asymmetric case Tl — > and 
T R — > 1 we shall bosonize only the electrons in the right 
contact. Then the unperturbed fixed-point Hamiltonian 
Hq + He is given by 



Ho = J2 £ ka{ a a ka + e P a : 



ka 



pa 



Hr; — Er 



-Vfe(o) 

IT ' 



N 



(46) 



(47) 



Here the operators au a and a pa correspond to the left- 
contact electrons in the lead and in the dot, respectively; 
the bosonic fields ^Ra an d n Ra model the electrons in 
the right contact. An integer- valued operator h describes 
the number of electrons transferred to the dot through 
the left barrier. 

The perturbation consists of two parts: weak tunnel- 
ing through the high barrier in the left contact, and weak 
scattering on the small barrier in the right contact, 



H' L = J2(v t al a a pa F + v* t al a a ka F r ), 



kpa 



H' R = -D\r R \y2coa[2<f> Ra (0)]. 



(48) 
(49) 



We have chosen not to write the operator n of the number 
of electrons transferred into the dot through the barrier 
in terms of the fermionic operators as n = ^ CLp a a pcn but 
to treat it as an independent variable. Thus the tunnel- 
ing Hamiltonian (|4^) acquires the charge-lowering and 
raising operators F and F^ defined as 



[F, n] = F. 



(50) 



One possible representatior£3 of these operators is given 



by F 



and h 



■ dx ' 



where x has the meaning of a 



phase operator conjugated to n. 

Our goal in this section is to find the conductance of 
the system, assuming that the perturbations H' L and H' R 
are small. At the first step we will utilize the smallness 
of the amplitude v t of tunneling through the left barrier, 
and find the expression for the current in the second order 
in v t . The current can be defined as the time derivative of 
the number of electrons that have tunneled through the 
barrier, 7 = —e^r — jr[n,H}. Using the commutation 
relation (p30|), we get 



I = i^Yj [y* t al a a ka F ] - v t a\ a a pa F 



kpa. 



(51) 



To find the conductance through the dot, we will sub- 
stitute the above current operator into the Kubo for- 
mula (E3). This requires the calculation of the correla- 
tor ([/(t), 7(0)]). Since 7 cx \v t \, up to the second order 



in tunneling matrix elements we can average [7(f), 7(0)] 
over the equilibrium thermal distribution of the system 
without tunneling, H' L = 0. Then the averaging over the 
fermionic degrees of freedom is straightforward, and we 
get 



G = -iG, 



ttT 2 



tK(t)dt 



Ti 2 J-oo sinh 2 [irT(t - iS) /%] 



(52) 



Here we have introduced the conductance of the left tun- 
nel barrier, Gl = {^e 2 /K) Y^k P \ v t\ 2 ^i. € k)8{e p ), and the 
correlator 



K(t) = (F(t)F^0)) = (Ft(t)F(0)). 



(53) 



The equality of the two correlators in Eq. (53) follows 
from the fact that the symmetry transformation n — ► — n, 
4>Ru — * — ^Raj&nd N — > —N does not affect the Hamil- 
tonian @, @, and @, but changes F+ <-> F. 

It is instructive to apply the formula (^2|) to the non- 
interacting case, Ec = 0. In this limit the operators F 
and F^ commute with the Hamiltonian (^6|) and (|4S|). 
Thus the correlator (^) is unity, and we get G = Gl- 
This is the expected result, since we assumed Gl 0, 
and the conductance of the system must be determined 
by that of the weakest link. 

Let us now find the conductance at non-zero charging 
energy Ec, but in the absence of the scattering in the 
right contact, r R = 0. In this case the correlator (|5^ ) 
is no longer trivial, because the operators F and F^ do 
not commute with the interaction term (|47|). When an 
electron tunnels into the dot, n — > n + 1, the system is 
shifted from its equilibrium state. To accommodate the 
additional electron brought by the tunneling process and 
to return the system to the equilibrium state, one elec- 
tron must leave the dot through the right contact. Thus 
the time evolution of the operators F and F> becomes 
non-trivial. 

To find the correlator K(t), we will perform the fol- 
lowing unitary transformation of the Hamiltonian: 



U = e™ e , 
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[Tl R1 (y)+U Rl (y)}dy. 



(54) 
(55) 



This transformation shifts the bosonic fields (f> Ra and 
adds a phase factor to the charge-lowering operator F, 



^Ra(x) -> (f>Ra(x) ~ - fl, 



F -> Fe 



/<-) 



(56) 
(57) 



One can easily see that the kinetic-energy term ( f46| ) 
is invariant under the transformation (|54|), whereas the 
charging term ( f47| ) transforms to a n-independent form, 



77 r 



U^H C U = E c 



<f> Ra {0)-N 



(58) 
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Upon this transformation, the operators F and F^ com- 
mute with the Hamiltonian (|4^) and but they ac- 
quire non-trivial phase factors, see Eq. (p7|). This enables 
us to find K(t) as the correlator of the phase factors, 
K(t) = K e (t) ee ( e i®(t) e -ie(o)). Since the phase 6 is 
linear in bosonic variables, the calculation of the correla- 
tor Kq is straightforward: 



#e(t) = exp{-< 

7T 2 T 



[e(o)-e(t)]e(o))} 
1 



2i~fE c sinh[7rT(i - iS) /%] ' 



(59) 



The last equality is valid asymptotically at low energies, 
T,h/t <C Ec- The substitution of K — Kq into the ex- 
pression for the conductance (|52| ) leads to a linear tem- 
perature dependence of the conductance: 



G — Gl 



8~fE c 



(60) 



Physically the origin of the linear suppression of the 
conductance atJow temperatures ( |60| ) is the orthogonal- 
ity catastrophe!^ Unlike the non-interacting case Ec = 
discussed above, the charge fluctuations in the dot are 
now suppressed at energies below E C - Effectively one can 
interpret the effect of the charging energy as a hard-wall 
boundary condition for the wavefunctions of the right- 
contact electrons with energies within the band of width 
W ~ Ec near the Fermi level. When an electron tunnels 
through the left barrier, the system must lower its charg- 
ing energy by moving one electron through the right con- 
tact. Since the two spin channels are completely symmet- 
ric, each mode transfers the charge q = I . According to 
the Friedel sum rule, this leads to a change in the bound- 
ary condition corresponding to an additional scattering 
phase shift 6 — ±^2 = ±-|. A sudden change of the 
boundary conditions creates a state with a large numbpc 
of electron-hole excitations, which is nearly orthogonalE3 
to the ground state of the system. This orthogonality 
leads to a power-law suppression of the tumipliag den- 
sity of states v(e) cx e x , with the exponentEJ~E3 deter- 
mined by the phase shifts in all the electronic channels, 
X = XX'V 71 ") 2 - I n our system, each spin mode gives two 
channels — in the dot and in the right lead — leading to the 
total of 4 channels. Thus the exponent x is unity, and 
the density of states v[e) cx e. The linear suppression 
of the tunneling density of states results in the T-lincar 
conductance (p0|). 

The linear temperature dependence of conductance 
( |60| ) appears to contradict to the T 2 law for the inelastic 
co-tunneling, Eqs. (||), ( pi] ) and (j39|). The reason is the 
absence of the barrier in the right junction assumed in 
the derivation of Eq. ([30|). We will now show that the 
presence of arbitrarily weak barrier in the right contact is 
crucial (see also Appendix A3) and leads to the quadratic 
temperature dependence of the conductance at T — ► 0. 

The weak backscattering in the right contact is de- 
scribed by the term (E3) in the Hamiltonian. As it follows 



from Eq. ([56]), under the unitary transformation (Q) the 
scattering term ( ff9|) acquires additional sign: 



H' R = {-lf-D\r R \Y d ^<t> 

TT ' * 



a. 



(61) 



The operators F and F' no longer commute with the 
Hamiltonian: they commute with H + He and anti- 
commute with H' R . Therefore the correlator ([53]) is no 
longer equal to Kq, but rather K(t) = K&(t)Kp(t), 
where 



K F (t) = (F(t)F\0)) 



(62) 



Here the averaging (. . .) is over the equilibrium ther- 
mal distribution of the transformed Hamiltonian, H = 
H + He + H' R . 

In the absence of the backscattering in the right con- 
tact, Kq is given by Eq. (|59"1), and K F (t) = 1. When 
a weak backscattering is added, both correlators are af- 
fected. It is clear, however, that Kq is not modified sig- 
nificantly. Indeed, it follows from the above mentioned 
analogy with the orthogonality catastrophe that the time 
dependence (59) is determined solely by the charge trans- 



ferred through the right junction, which in our case is al- 
ways e — the charge of the electron brought into the dot 
through the left barrier — and is not affected by the pres- 
ence of the right barrier. Therefore the only change in 
KQ(t) caused by the weak barrier can be a small modi- 
fication of the prefactor in Eq. ( |59]) . On the other hand, 
we will show that K R (t) is modified significantly, which 
leads to the quadratic temperature dependence of the 
conductance. 

To find Kp(t) we first note that the Hamiltonian 
H = H + He + H' R is identical to the one considered 
in Sec. HI and in Ref. |lq. The difference with the prob- 
lem considered in Sec. |III| is that instead of two modes 
corresponding to two different contacts we now have two 
spin modes in the same (right) contact. Thus the two 
scattering amplitudes are equal and given by |t\r|. As we 
have seen in Sec. HI, the presence of the backscattering 
shows up at low-energy scales e ~ \r R \ 2 Ec <C Ec, where 
one can effectively integrate out the mode related to the 
charge fluctuations in the dot, 4>R\ + 4>Rl- The result- 
ing Hamiltonian is expressed in terms of the spin field 
0.Rs = (0at — < ^ ) i?4)/v / 2 and has the form: 

H= ^¥j {-[WflsW] 2 + rflL(z) J dx 

V^0iu(O)l. (63) 



+ (-!)' 



8-fE c D 



r R | cos 7r N cos 



We can now fermionize the Hamiltonian using a represen- 
tation of the operators e ±4v/2 ^ Hs in terms of the fermion 
creation and annihilation operators in a way identical to 
the one leading from the Hamiltonian (|2^) and (|3(]) to 
Eq. (^). The final Hamiltonian is again quadratic in 
fermion operators, but also depends on the operator n 
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of the number of particles transferred through the left 
barrier, 

" 2 



+ (-l) fl -yJ^hv F E C \tr\ COS nN 
IT 

x(c + c<)[tj;(0)-^(0)}. 



(64) 



Here we presented a Majorana fermion operator similar 
to 77+ in Eq. ( |35| ) as a sum of fermion creation and an- 
nihilation operators, rj = c + . This representation is 
more convenient, because we can now replace the oper- 
ators F and in the definition (jf^) of the correlator 
Kp by F s = F\ = \ — 2c 1 * c. Indeed, these new oper- 
ators commute with the first term in the Hamiltonian 
(|64|), anti-commute with the second one, and possess the 
property F a FJ ^^j—FF 1 * = 1. Therefore one can find the 
correlator Kp ast 



K F {t) = ([1 - 2ct(t)c(t)][l - 2c t (0)c(0)]>. 



(65) 



The actual calculation of the correlator (|65]) is now 
straightforward, as the Hamiltonian (|64|) is quadratic in 
the fermion operators. The result for the correlator Kp 
has the form 



K F {t) = 



2r, 



AEt 



dE 



(E 2 



Here we have introduced the low energy scale 
T R = ^E c \r R \ 2 cos 2 nN, 



(66) 



(67) 



originating from the presence of a weak scatterer in the 
right contact. In the weak-scattering limit, Tr — ► 0, 
Eq. ( |66| ) reproduces the result Kp(t) = 1 mentioned 
above. However, in the opposite limit T, %/t <C Tr, the 
correlator Kp has the same non-trivial time dependence 
as the correlator K@ (t) given by Eq. (p39|) , 



K F (t) 



2T 



1 



iF R sinh[7rT(t - iS)/h] ' 



This additional time-dependent contribution to K(t) 
leads to the quadratic temperature dependence of the 
linear conductance at T -C Tr. 

To find the conductance at temperatures T ~ Tr, we 
can now substitute the correlator K(t) = K®{t)Kp{t) 
into Eq. (|52l), which leads to the following expression: 



G = 



{E/Tf 



TG, 



8jE c E 2 + T 2 R cosh 2 (E/2T) 



dE. 



(68) 



Equation (68) is the central result of this section. It 
is valid for any T/Tr, if both temperature and Tr are 
much smaller than the charging energy Eq- The condi- 
tion Tr <C Ec means |7"r| <C 1, i.e., the scattering in the 
right contact must indeed be weak. 
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FIG. 3. Conductance (^8|) as a function of the dimension- 
less gate voltage iV for the asymmetric case, \rn\ = 0.4. The 
three curves are calculated for Ec /T = 10 (a), 40 (b), and 
160 (c). 



In the limit T ^> Tr, Eq. ( pq ) reproduces the pre- 
viously obtained linear temperature dependence ( |60| ) in 
the absence of the scattering in the right contact. It is 
important to note, however, that the linear temperature 
dependence is obtained not only at tr = 0, but also in 
the presence of the barrier if N is half- integer, i.e., at the 
centers of conductance peaks. 

Away from the centers of the peaks, the conductance 
has a stronger temperature dependence, G oc T 2 . Indeed, 
at T -C Tr the expression (|68|) reduces to 



G 



2tt 2 T 2 
37 E C T R 



G, 



(69) 



Thus the quadratic temperature dependence of the lin- 
ear conductance is restored for the off-peak values of the 
gate voltage. The evolution of the conductance peaks is 
illustrated in Fig. |[ 



V. SCALING APPROACH TO INELASTIC 
CO-TUNNELING 

In the sections ||-4Tv| we have considered several cases 
when the inelastic co-tunneling can be studied analyti- 
cally even though one or both contacts are in the strong- 
tunneling regime. A common property of all these re- 
sults is that at T — > the conductance behaves as T 2 
away from the resonance values of the gate voltage, i.e., 
at N ^ ±i,±|,... The same temperature dependence 
was obtained earlierEI for the off-resonance conductance 
in the weak-tunneling case, see Eq. (0). In this section we 
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study the general scaling properties of our system and ar- 
gue that the quadratic temperature dependence of linear 
conductance is universal and should hold for arbitrary 
barriers. 

We will describe our system by the bosonizcd Hamil- 
tonian (|19|)- (pi]) in the spinless case and (^o|) -(^2|) in the 
presence of spins. In the absence of the backscattering 
the Hamiltonian is quadratic, i.e., \r^\ = \r R \ = is 
a fixed point of our Hamiltonian. In analogy with the 
Kondo problem we will call it the strong-coupling fixed 
point. To investigate the stability of this fixed point, one 
should find the scaling dimension of the backscattering 
term or @. 

We shall start with the spinless case, where at T = 
and t 3> h/Ec a simple calculation gives 



nr (N) 



2nit 

-3 



W{\r L \ 2 + \r R \ 2 - 2\r L \\r R \ cos2^V) 
4iriE c t 3 



(H'(t)H'(0)) = 



with Tq(N) defined by Eq. (|26|). The long-time asymp- 
totics of the correlator (H'(t)H' (0)) is given by the first 
term, which is proportional to t~ 1 , unless To = 0. Thus 
the dimension of the perturbation H' is \ < 1, and the 
strong-coupling fixed point is unstable. At a half-integer 
N and in a perfectly symmetric case, \tl\ = \r R \, the 
parameter Tq vanishes, and the dimension of the opera- 
tor H' becomes | > 1. In this case the strong-coupling 
fixed point is stable. In the language of the correspond- 
ing 2-channel Kondo problem this means that the strong- 
coupling fixed point is stable only in the absence of mag- 
netic field and channel anisotropy. The leading irrelevant 



operator with dimension ~ gives rise to the weak temper- 
ature dependence of the peak heights in the symmetric 



spinless case, see Appendix A 1 



A similar analysis can be performed in the case of elec- 
trons with spins, described by the Hamiltonian (fiol)-([42"l). 
The correlator of the perturbation operators now has the 
form 



(H'(t)H\0)) = (\r L \ 2 + \r R \ 



c 



3/2 



meaning that H' has dimension | < 1 at any N and any 
barrier strengths. Thus as expected the strong-coupling 
fixed point in the 4-channel case is unstable. At a half- 
integeiJV this statement has a well-known analog in the 
theorycZI of the multichannel Kondo problem, where the 
strong-coupling fixed point is unstable even in the ab- 
sence of the magnetic field. 

The fact that the strong-coupling fixed point is unsta- 
ble means that as the temperature is lowered, the weak 
backscattering in the contacts becomes stronger. It is 
natural to assume that the effective barriers grow in- 
definitely, and the system approaches the weak-coupling 
fixed point. To test this hypothesis we will consider a 
system where the two contacts are in the weak-tunneling 



regime and study the stability of the corresponding fixed 
point. 

As we saw in Sec. || and IV, in the weak-tunneling 
case it is convenient to describe the system in terms of 
the original fc rmion ic variables. Generalizing the Hamil- 
tonian of Sec. IV B to the case of both barriers being in 
the weak-tunneling regime, we introduce the Hamiltonian 
H = H + H C + H' L + H' R , where 



ka pa 

y ] e q a qa a qa + ^ ] £ s tt sa tt si 



qa 



H c =E c (n L +n R -N) 2 , 

H L = J2( v L a la a p a F L + v L a la a ka F l)i 



kpa 



H R = J2( v R a la a qa F R + v R a \a a sa F R )- 



(70) 

(71) 
(72) 

(73) 



qsa 



Here ak a and a sa describe the electrons in the left and 
right leads, respectively; a pa and a qa correspond to the 
electrons in the dot coupled to the l eft an d right contacts. 
Similarly to the formalism of Sec. IV B, we have intro- 
duced the operators n l and n R of the number of electrons 
tunneled through the left and right contacts and the cor- 
responding raising and lowering operators F L , F^, F R , 
and F R in the tunnel Hamiltonians (72) and (|73|). 

Away from the resonance values of the gate voltage, 
N 7^ ±i ±1,..., the states with different nj, and n R 
have different energies because of the Coulomb term ( |7l|) . 
Since the tunnel Hamiltonian given by ([72]) and (|73|) 
changes n;, or n R , it does not describe the transitions 
between the low-energy states of the Hamiltonian. Thus 
when considering the low-energy properties of the sys- 
tem, one should use the tunneling terms @ and (^§ to 
construct the perturbation which is not accompanied by 
the change of the total charge of the dot. In the simplest 
case, this can be done by selecting the terms containing 
equal number of operators Fi and Fj in the second-order 
perturbation 



H! 2 = (H' L + H' R ) 



1 



En — Hn — H t 



c 



-{H' L + H' R ). (74) 



Here the denominator accounts for the energy of the vir- 
tual state created by adding (or removing) an electron to 
the dot. At very low energies the only important contri- 
bution to the denominator is the change of the electro- 
static energy of the system. 

At iV in the interval (—5,5) a typical representative 
of the variety of terms contained in Eq. (^41) has the form 



A = -v* L v R 



E C {1-2N) E c (l + 2N) 



xF F. 



• j ± R a pa a ka J! ' 
kpa qsa' 



(75) 



12 



The terms included in Eq. ( |75| ) correspond to the pro- 
cesses where one electron tunnels into the dot through 
the left junction and another one escapes from the dot 
through the right junction. Depending on which of the 
two tunneling events is executed first, we get different 
energies of the virtual states, corresponding to the two 
denominators in Eq. (f75|). 

To find out whether or not the weak-coupling fixed 
point is stable, one should evaluate the scaling dimension 
of operator A and other similar contributions in Eq. (74). 
A straightforward calculation gives 



(A(t)AH0)) = 



(27re 2 ) 2 



1 



-1 2 



Ec(l 
ttT 



+ 



2N) E c (l + 2N) 



swh[jrT(t-i8)/h} 



(76) 



At T — ► the correlator decays as 1/i 4 , and therefore the 
dimension of the operator A is 2. One can also check that 
the other contributions to Eq. (174) have the same dimen- 
sion, whereas the higher-order terms have even higher di- 
mensions. Thus the weak-coupling fixed point is stable. 

So far we have shown that away from the reso- 
nance values of the gate voltage, N ^ ±i,±|,..., the 
strong-coupling fixed point is unstable whereas the weak- 
coupling fixed point is stable. We now conjecture that 
there are no fixed points at intermediate coupling. Then 
at low energies the system always evolves towards the sta- 
ble (weak-coupling) fixed point. One should also expect 
that near this fixed point the temperature dependence of 
conductance is determined by the dimension of the lead- 
ing irrelevant perturbation, which is a universal charac- 
teristic of the problem. Indeed, one can easily show that 
the leading contribution to the current operator I = en^ 
is given by 



I 



n 



(A -At). 



(77) 



We now substitute this current operator into the Kubo 
formula (|2^) , and use (|76| ) to find the following expression 
for the conductance: 



G 



■kHGlGr 



1 



1 



-1 2 



E c {l-2N) E C {1 + 2N) 



(78) 



This expression for the conductance in the weak- 
tunneling limit coincides with the well-known result for 
the inelastic co-tunneling,B and reproduces Eq. (||) near 
the resonances, \ ± N <C 1. 

It is important to note that in the derivation of Eq. ( |7g| ) 
we neglected all the contributions to the current opera- 
tor with dimensions d > 2. Such terms would lead to the 
additional terms in the correlator ([76|), which are propor- 
tional to {ttT/ smh[nT(t-iS)/h]} m with m > 4. Clearly, 
the corresponding corrections to the conductance would 
behave as SG cx T m ~ 2 < T 2 . Thus at T -> the tem- 
perature dependence of the conductance is given by the 



leading irrelevant perturbation. As we have seen, the uni- 
versal scaling dimension d = 2 of this perturbation trans- 
lates into the universal temperature dependence of the 
conductance G cx T 2 . Therefore this temperature depen- 
dence is specific not only for the cases of weak or strong 
tunneling, but should hold for any barrier strengths. 

The above treatment of the stability of the weak- 
coupling fixed point was limited to the off-resonance val- 
ues of the gate voltage. We will now demonstrate that the 
properties of the system are completely different at half- 
integer N. Consider, for example, the case of N = \. 
Then the states with the dot charge equal to and e 
have the same energy. As a result the terms in the tun- 
nel Hamiltonians (|7|) and ( |73| ) responsible for the tran- 
sitions between these two states no longer lead to the 
large increase in the energy of the system and should 
not be discounted. It is easy to check that the dimen- 
sion of such operators is d — 1, i.e., the tunnel Hamilto- 
nian is a marginal operator. This conclusion is actually 
obvious, because the same calculation of scaling dimen- 
sion must be valid in the absence of the charging effects, 
Ec = 0, where we do not expect the tunneling ampli- 
tudes to be renormalized as the temperature is lowered. 
However, the simple calculation of the scaling dimension 
is not sufficient to determine the stability of the weak- 
coupling fixed point. To make such a determination, one 
has to study the scaling equations for the tunneling am- 
plitudes including the higher-order terms in vl and vr. 
We have performed this analysis in Sec. where we saw 
from the mapping to the multichannel Kondo problem 
that the tunneling amplitude is in fact a marginally rel- 
evant perturbation, i.e., it grows logarithmically at low 
temperatures. Therefore on resonance the weak-coupling 
fixed point is unstable. 

The on-resonance scaling properties of our system are 
strongly affected by the presence of electron spins and by 
the symmetry of the two barriers. Let us first consider 
the case of symmetric barriers. In the spinless case we 
saw that on resonance the strong-coupling fixed point is 
stable if the barriers are symmetric. Thus one should ex- 
pect the system to always renormalize towards this fixed 
point. Then the conductance must be on the order of 
e 2 /h. This is confirmed by the explicit calculation in the 
strong-tunneling case, Eq. (p8[). 

In the presence of spins, both weak- and strong- 
coupling fixed points are unstable on resonance. Thus 
one should expect that there must be a stable fixed point 
at some intermediate coupling strength. This hypothe- 
sis is also confirmed by the analogy with the 4-channcl 
Kondo problem discussed in Sec. ||. Indeed, it is well- 
known that the low-energy properties of the Kondo prob- 
lem with the number of channels exceeding 2 is governed 
by an intermediate-coupling fixed point, which gives rise 
to a specific non-Fermi^auieLbehavior of the suscepti- 
bility and specific heat.E3cJ~E3 In Appendix |b] we show 
how one can study the properties of the intermediate- 
coupling fixed point starting from the vicinity of the 
strong-coupling point. We discover that on resonance 
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the thermodynamics of our problem is indeed identical 
to that of the 4-channel Kondo problem. 

The technique presented in Appendix [b|, however, does 
not enable us to find the low-temperature conductance G. 
Nevertheless one can estimate the order of magnitude of 
G at the intermediate-coupling fixed point in the follow- 
ing way. From the fact that the strong-coupling fixed 
point is unstable and from the result (|4|) of the pertur- 
bation theory near this point we know that the conduc- 
tance must be smaller than e 2 /2nh. On the other hand, 
at G <C e 2 /h the weak-coupling treatment of Sec. || is 
applicable and predicts the growth of conductance at 
T — > 0. Therefore the limiting value of conductance 
must be of the order, but smaller than e 2 /2nh. We fur- 
ther conjecture that the peak conductance is a universal 
characteristic of the intermediate-coupling fixed point, 
see Appendix |b], meaning that in the limit T — > the 
peak conductance does not depend on the heights of the 
barriers. 



The scaling analysis in the spirit of Ref. 23 enables us 
to make some conclusions about the shape of the peaks 
in G(N). In the spinless case we see from Eq. ( |3^ ) that 



at the tails of the peak at N 
i\ 4 



the conductance de- 



cays as G cx 1/ (N — 5)^- This result can be interpreted 
within the scaling picture as follows. A small deviation 
of N from ^ is identified as a magnetic field h in the 
appropriate multichannel Kondo problem, see Eq. (0). 
Therefore such a deviation is a relevant perturbation with 
scaling dimension d = 2/(2+ k), where k is the number of 
channels .lEI In the spinless case we have k = 2 and d = i. 
This means that as the temperature is lowered, the mag- 
netic field grows as h(T) oc (N - i)T -1 / 2 . We now 
conjecture that near the fixed point, at T — > 0, the con- 
ductance is a universal function of h only, G = G(h(T)). 
As we discussed above, away from the center of the 
peak G cx T 2 at T — > 0. Therefore the magnetic 
field dependence of the conductance must be given by 
G cx l/h 4 cx T 2 /(N - i) 4 , in agreement with Eq. Q. 



We can now apply the same argument to the case of elec- 
trons with spins, where the number of channels k = 4 and 
the scaling dimension of the magnetic field is d — -| . The 
temperature dependence of the renormalized magnetic 
field is consequently h(T) oc (JV - \)T~ 2 ^. Therefore 
near the intermediate-coupling fixed point the condition 
G oc T 2 at T -► demands G oc h~ 3 oc T 2 /(N - ±) 3 . 
Thus we expect the tails of the conductance peaks to 
decay as l/(iV — ^) 3 . 

So far we discussed the case of symmetric barriers. The 
on-resonance behavior of the system is very different if 
the two barriers are not identical. This can be readily 
seem from the analogy to the Kondo prq blpr n. From the 
theory of the multichannel Kondo modelEIrEZI it is known 
that the asymmetry between the channels is a relevant 
perturbation. In terms of the Coulomb blockade prob- 
lem this means that on resonance a small asymmetry 
will grow and eventually the smaller transmission coeffi- 
cient, say Tl, will scale to zero, whereas the larger one, 



Tr, must approach the fixed point of the problem with 
half the number of channels of the original one. Thus for 
both 2- and 4-channel cases Tr — > 1. For the most re- 
alistic case of electrons with spins, near this asymmetric 
fixed p oint t he system is identical to the one considered 
111 Sec. |IVB[ Thus we conclude that for asymmetric sys- 
tems (i) the on-resonance conductance vanishes linearly 
in T, and (ii) the shape of the resonances in G(N) is 
universal and given by Eq. (|6§|). 

The on-resonance scaling properties of our problem are 
summarized in Fig. ^|. 




(b) 




r R 



FIG. 4. Scaling diagram of the two-barrier system in 
the cases of (a) spinless electrons and (b) electrons with 
spins. The gate voltage corresponds to a resonance value, 
N = ±|, ±|, .... In the symmetric case, TL = Tr, the system 
scales to either strong- or intermediate-coupling fixed point, 
depending on the number of channels. In an asymmetric case 
the system approaches the limit where the weaker barrier dis- 
appears, while the stronger one becomes very high. 
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VI. CONCLUSION 

In this paper we have studied the transport through a 
quantum dot connected to two leads by quantum point 
contacts, Fig. [I]. -Unlike in conventional theories of 
Coulomb blockade,BB we concentrate on the regime of 
strong tunneling, where one or both contacts are close 
to perfect transmission. We find the linear conductance 
in a wide interval of temperatures, limited by the level 
spacing in the dot from below and by the charging energy 
Ec from above. 

We have shown that as long as there is any backscat- 
tering in the contacts, i.e., both transmission coefficients 
Tl.r < 1, the Coulomb blockade peaks in the conduc- 
tance as a function of the gate voltage should be ob- 
served. Between the Coulomb blockade peaks the tem- 
perature dependence of the conductance at T — > is de- 
termined by the weak-coupling fixed point. As a result we 
get the T 2 dependence of ihe conductance, which coin- 
cides with the known resulirJ for the inelastic co-tunneling 
mechanism. 

Our results show that the behavior of the conductance 
near the peaks depends qualitatively on the presence of 
electron spins. In the spinless case the strong-tunneling 
fixed point is stable if the barriers are symmetric, and 



thus the technique developed in this paper allows a com- 
plete solution of the problem. The resulting peaks in con- 
ductance have the height e 2 /AttH, with the shape given 
by Eq. (|38l). In the presence of the spins, the strong- 
tunneling fixed point is no longer stable, and our tech- 
nique does not allow a complete solution of the prob- 
lem. Using a useful analogy with the multichannel Kondo 
problem, we have concluded that the peak conductance 
must still be of the order of e 2 /%. In addition, we stud- 
ied the peaks in conductance in the case of asymmetric 
barriers. The peak conductance in this case is linear in 
temperature at T — ► 0, and its shape is given by Eq. (|6g|). 

The strong-tunneling regime of Coulomb blockade wag 
explored to some extent in several recent experiments.LTtl 
Our results, are in qualitative agreement with the 
experiment B We believe that our predictions can be eas- 
ily tested in similar experiments. 
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APPENDIX A: PERTURBATION THEORY FOR WEAK REFLECTION 



In this appendix we derive the leading-order correction to the conductance due to the weak reflection at the quantum 
point contacts. We shall use a path integral method to calculate the current-current correlation function at imaginary 
frequencies and then perform analytic continuation to the real frequencies. 



1. Spinless fermions 

In the spinless case our problem is similar to the resonant tunneling problem studied in Ref. ^l|. Thus we can apply 
the method developed in Refs. |2l| and [38] to calculate the conductance perturbatively with respect to the reflection 
coefficients. 

The current operator is given by I — (e/ \f2n)dt4>i{Q, t). Then the Kubo formula ( p2| ) for the conductance can be 
rewritten as 

e 2 T 

G = - „ lim u> lira (<fii(-iuj n )<fri(iujn)), (Al) 

2-!T 2 in iuj n -*uj+i0+ 

where u) n — 2imT/h and 

f h/T 

M*o»n) = / e^"> / (0,r)dr. (A2) 
Jo 

The thermal average $/(«o; rl ) = (4>i(~iuj n )(j)i(iuj n )) is calculated as 

1 S 2 Z W 



Zy) dJ(lLO n )dJ{-lUJn) 



(A3) 



where the generating functional is given by 
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St 



— - + J{—iu) n )4>i{iu n ) 



Si 



2vr Jo 



h/T 



dr / dx 



— {dr^f + V F (d x ^) 2 + —(d T ^c) 2 + V F {d x( j>c) 2 



2E C 



h/T 



dr 



<MO,r) 



nN 

71 



D [ h ' T 



dr ( | r L | cos {V2 [<^/(0, r) - C (O, r)]}+ |r«| cos{V2 [0/(0, r) + O (O, r)]} 



(A4) 



Expanding exp(— S\/K) in powers of |r| and performing Gaussian integrals, we get 



2\cj n \T 2^2 

7T 4 7 T 2 



[|r L | 2 + M 2 + 2|r L ||r fl |cos(2^V) 



h/T 1 - cosw„r , 

rrrdr 



sin(7rTY/7i) 



[K| 2 + M 2 - 2|r L |Mcos(27r7V)] 



In this calculation we have used the following relations: 

T e -H\u>' n \/D i 
jX Ml (l-coso; n r) = -ln 



fi/T 



■ sm 



1 — cos W n T 
sin 3 (7rTT/?i) 

7tTt\ 



dr + 0(|r|' 



ttL> 



T e -?»Kl/g _ 1 

h *f \iv' n \ + (2Ec/Trh) ~ tt U \2~fE c 

7T 3 T 2 



h^-f \^ n \ + {2E c /irK) 



cos Cl)„ t 



4£ 2 sin 2 (ttTt/Tl) 



(A5) 

(A6a) 
(A6b) 
(A6c) 



where we have neglected higher-order terms in T/Ec- Note that Eqs. (A6a) and (A6c) are correct only for r, fiT 1 — 
r > h/D. 

In calcu lating the integrals over r in Eq. ( A5), we introduce a short-time cutoff r c ~ ?i/-D, which is omitted in 
Eqs. ( A6a ) and ( A6c ), and modify the contour of integration: 



h/T 



1 — COS U n T 

sin" '(wTt/H) 



dr = 



h/T-T 



1 - e* 1 



sin^TrTr/Ti) 



dr = 



OO ^ _ g UJ n {T c +tt) 

sin" [ttT(t c + it) /h] 



dt 



OO ^ _ e iLU n (-T c + it) 

q sin 1/ [7rT(r c — it)/h] 



By replacing iw„ by w and taking the limit uj — > 0, Eq. (A7) is reduced to —iujF(T, i>), where 



F(T,u) = 



it 



rdt = 



dx 



/5Fr(|) 



sia v [7rT(r c + it)/n] 2ttT 2 7„ oc cosh" x 2irT 2 T(f + |) 



dt. (A7) 



(A8) 



From Eqs. (Al). (A5), and (A8) we finally obtain 



G = ^j 1 - ^[Kl 2 + l^l 2 + 2M|rfl|cos(27ri\0] - ^[\r L \ 2 + \r R \ 2 - 2|r i ||r fl |cos(27rJV)] ). (A9) 



2ttT 



16E C 



The term proportional to T/Eq originates from the weak fluctuations of the charge of the dot at low frequencies. At 
a half-integer TV it is related to the leading irrelevant operator near the fixed point of the 2-channel Kondo problem. 
The scaling dimension of this operator is d = | (see Sec. |y|), yielding the temperature dependence SG oc T 2d ~ 2 = T. 



2. Electrons with spins 

The calculation proceeds in parallel to the spinless case. We first define the new sets of phase fields, 

<t>Lc = ^(^T + <f>I>i), ^Ls = -^=(0Lt - <f>Ll), 0Rc = "^(^T + <f>Rl)i <t>Rs = ^ (<^T ~ <i>Rl)> ( Al 0) 

where <pL C and 4>Rc describe the charge density fluctuations, and 4>Ls and cj>R S correspond to the spin density fluctua- 
tions. We further define their momentum conjugates IIl c , TLl s , ^Rc and Hr s in the same way. 
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The current operator in this basis is I — (e/\/2ir)dt[<f>Lc(0,t) + 4>R C (Q,t)], and accordingly the Kubo formula for 
the conductance is 



G 



e 2 T 



lim ui lim $ c (iu; n ), 



7J 11111 UJ 11111 



where $ c (iw n ) is given by 



5 2 zf 



SJ(~iuj n )SJ(iu! n ) 



.1=0 



-jr(S 2 + Sj) 



with the generating functional, 

Zj 2) = jT>(f>Lc(x,T) jT>(f)L S (x,T) Jt>(/) Rc (x,t) Jv<j) Rs (x,T)eyLp 

— {d T (j> Lc ) 2 + v F (d x <j> Lc ) 2 + —(d T <p L s) 2 + v F (d x <f> Ls ) 2 
v F v F 

+ — (dr^c) 2 + v F {d x <p Rc f + — {d T cj)B. s ) 2 + v F {d x <j> Rs ) 

v F 



(All) 



(A12) 



(A13) 



5, = 



h 

2Wo 



h/T 



dr / dx 



2E r 



r h/T 



dr 



v F 



+ — J rfr||r L |cos[V20 Lc (O,T)] cos[V20 Ls (O, t)] +|r Jt |cos[V20 fl , c (O,T)] cos[V20 fls (O, r)]|, (A14) 

Sj = h]Tj(-iU} n )[<f>Lc(iUn) + 0Rc(iUn)}- 



Up to the order \r\ 2 the correlation function Q c (ioj n ) is obtained as 



rh/T 

~-^VlTE^{\rL\ 2 + \rn\ 2 ) / - 



irk 

\u) n \T flUJ 2 



1 — COSW n T 



sin 3 / 2 (7rTr/^) 



■dr, 



(A15) 



where higher-order terms in T/Ec are neglected. Note that there is no term proportional to |rj,||rij| in Eq. ( |A15| ) 
because (cos</>l s (0)) = (cos4>r s (0)} = 0. This is an important difference between the two-channel case and four- 
channel case, implying that the strong-tunneling limit is unstable even on resonance in the latter case. From Eqs. (|A7| ), 
@, (|Al~l"l), and (|AT5|) we can readily obtain Eq. (El). 



3. Asymmetric limit 



We have discussed the case of asymmetric barriers in Sec. IV B, where the conductance is calculated in lowest order 
in v t but in all orders in \rn\ by effectively summing up all the most divergent terms in perturbation series. Here we 
shall calculate the leading term in the series, which diverges at low temperatures. 

We will calculate the conductance in a way slightly different from the one used in Sec. |IVB| . We first note that since 
the left barrier is very high, the left phase fields are pinned at the minima of the potential cos(\/20Lc) cos(\/2<^Ls) , 
i.e., (4>Lc: 0Ls) = (•\/27r(m+'|), y/2nnj or (^/2irm, \/27r(ri-\-h)), where m and n are integers. Thus the tunneling of 
an electron through the left barrier corresponds to a sudden change in the phase fields, say, (4>Lc, 4>l s ) = (tt/ \/2, 0) — > 
(v27r, 7r/v2)< The probability for this tunneling process is then given by 

P+ = J^! [°° ^mt/^-iHft/H^ieVt/H^ (A16) 
J — oo 

where v t is a constant proportional to vt, V is a voltage across the left barrier, and Hi and Hf are the bosonized 
Hamiltonian (|40|)-(p^) with the condition (4>Lc{0), 4>Ls(0)) = (tt/V^O) and (y/2n, 7t/a/2), respectively. The thermal 
average (. . .)j is taken for the Hamiltonian Hi. It is clear that the probability for the inverse tunneling process is 
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given by P_ = e eV / T P + . Since the current is I = 2e(P+ — P_) with the factor 2 coming from electron spins, the 
linear conductance is given by 



G 



h 2 T ' 



2\~ i 2 r°° 



h 2 T 



K (t)dt. 



(A17) 



It will soon become clear that the correlator Ko(t) is related to K(t) in Eq. ( p2[ ) as Ko{t) — {ttT/iD sinh[7rT(i — 
iS)/h]} 2 K(t). For the imaginary time t — it we can calculate the correlator using the path integral. To get the 
conductance we then need to perform analytic continuation, r — > it + 5, and calculate the integral ( ]A17[ ). 

The partition function of the system Zj | j = o is given by Eq ( |A13| ) . Since both the charging energy and the cosine 
potential depend on the phase fields at the point x — only, we can integrate outES the phase fields off the point to 
get the following effective action: 

ScS = — ^ [l < / , £c(«W„)| 2 + \<j)Ls(iu n )\ 2 + \4>Rc{iUn)\ 2 + 1 4>Rs (w n ) 1 2 ] 



+ - 



+ - 



2E C 

7T 2 

2D 



h/T 



dr 



4>r c {t) - <Plc{t) - 



h/T 



dr | \r L | cos [V2<f> Lc (rj\ cos[V20 Ls (t)] + \r R \ cos[V2<?W(t)] cos [V2<p Rs (t)] | . 
With this effective action we can readily calculate Ko(—ir) as 



Koi-ir) 



JV(j) Rc J V!p Rs exp{-j-S cS [>^j , , 4> Rc , cj> Rs ] } 

f T>(j) Rc J V<l) Rs exp{ - i S eS [4^1, <t>i/ a , 4>Ro , 4>Rs] } 



where 



# c V),0i°V)) = (t/V2,0) < r' < ft/T, 



(«&V),*kV))- 



27r,7r/V2) 0<t'<t, 
tt/\/2,0) t<t'<H/T. 



(A18) 

(A19) 

(A20a) 
(A20b) 



We shall first consider the non-interacting case, Eq = 0, to find the coefficient in Eq. ( A17 ). In this case we may 
ignore the fields <p Rc and 4> Rs . Thus the correlator is given by 



K (-it) 



exp{-±S e 4$l $1,0,0}} 



TFT 



Dsm^Tr/n) 



(A21) 



where we have used Eq. (A6a). Since we know that in this case the conductance is G = Gl, we rewrite Eq. (A17) as 



G 



G L D 



2 roo 



K (t)dt. 



2ttTiT _ 

With nonzero charging energy, Ko(—ir) is calculated up to the order |?\r| 2 as 



K (-ir) = 



7T 4 T 3 



1 



1 E C D 2 sin 3 (ttTt/K) 



1 



T R T 



dr\ 



h/T 



dr- 



2 — 



sin[7rT(T 2 — Ti)/K] J 



(A22) 



(A23) 



where T R is defined in Eq. ([37|). In this calculation we have neglected terms which are smaller in T/Ec such as 
TJ2^ sin(o; n Tj)[|a; n | + (2Ec/^h)]~ 1 . For < r < ^ the double integral ^(t) is calculated as 



pr T nh/2T 

W)=2 / . /M dr -2r / 
J sm(nTT /h) J T 



dr 



sm{irTT /h) T \2T 



ih , 

— — -t +2 



sin[7rT(r + it )/h] 



dt . (A24) 



To get the conductance we need the following integral: 
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m dt = 2 (jl 

sin 3 (i7rTt/n) \nT 



dx I dy 
'o 



cosh x cosn(x + y) \7ri , 



From Eqs. ( |A22[ ), ( |A23[) , and (|A25|) we get 



G = Gj 



TT T 

'WEc 



l-g[7C(3)-2] 



(A25) 



(A26) 



The leading correction due to the weak reflection in the right contac t div ergesEl at low temperatures unless N is 
half-integer. We note that for T > T R Eq. @) in fact reduces to Eq. ( |A26| ). 

The above result is for the case of electrons with spins. We can, of course, perform the same kind of calculation for 
the spinless case as well. Here we give only the final result: 



2tt 4 T 2 

G = °l^T^ ~ 47M cos(27riV)]. 



(A27) 



There are two important differe nces between Eqs. ( A26 ) and ( A27 ). First, even without the reflection in the right 
contact, the conductance (A27) is already proportional to T 2 . Second, the leading correction gives no additional 
temperature dependence. Thus the perturbation expansion in powers of |t\r| is well behaved for the spinless case, in 
striking contrast with the case with_spins. A similar difference between the two cases was found for the charge of a 
quantum dot with a single contact Ji3 



APPENDIX B: 4-CHANNEL KONDO FIXED POINT APPROACHED FROM THE STRONG 

COUPLING LIMIT 



In this appendix we demonstrate that, when N is half-integer (on resonance), our ID model is closely related to the 
4-channel Kondo model by tcansforming our bosonized Hamiltonian to the one that recently appeared in the study 
of 4-channel Kondo problem.L3 

First we rewrite the Hamiltonian (f40|)-(|42"|) with the fields introduced in Eq. ( A10| ) . We then introduce new fields 
tp Lc , <Pl s , Vrc, an d <Pfi s defined in the same way as ip + in Eq. each describing the even modes of 4>^ c , §l s , <p Rc , 
and 4> Rs . The Hamiltonian is now given by 



H 



VF 
4-7T 

E c 



{[V^ Lc (x)] 2 + [V^ Ls {x)f + [V^ Rc {x)f + [V^ s (x)] 2 } dx, 



Hv = ^Mcos[^ c (0)]cos[^ s (0)] + ^^cos[^ c (0)]cos[^(0)]. 



(Bl) 
(B2) 
(B3) 



With these new fields the current operator is written as 

evF 



2vr 



x=0 



Next we fermionize the bosonic fields as in (M) to get 



Hq = ihvp 



Hc=E c \J sgn(x) [i/>l e (x)i/> Le (x) - ^ Rc {x)i: Rc {x)\ dx - N' 
H V = -hv F \r L \ VLcVLs ty Lc (Q) - 4 C (Q)] [^(0) - 4.(0) 

-KvF\r R \v Rc V Rs [^(0) - 4j0)] [^(0) - 4,(0) 
/ = -ev F : 4 c (0)^ £c (0) + 4(0)^(0) : , 



dx, 



(B4) 

(B5) 
(B6) 

(B7) 
(B8) 
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where the rj's are Majorana fermions, the ifj's chiral fermions, and N' = N + const. Since the products of the Majorana 
fermions, rj Lc rj Ls and T) Rc r) R8 , commute with each other, and (VlcVls) 2 = (VrcVrs) 2 = ~ lj we can replace them by 
—i. We then introduce another set of fermions: 



i>2c( x ) = "4? \4 Lc ( X ) + #i?cO) 



i>s( x ) = 2 ^Ls{ x ) + 4 Ls ( x ) + + ilpR s ( x ) 

In terms of these fermions the Hamiltonian and the current operator are written as 



Hq = iHvf 



dx. 



H c = E c \ / sgn(x) U\ c (x)^ lc (x) - 4 2c ( x )^2c( x ) 



dx-N' 



H v = 



—%v F \r L 



I = —evF 



^ic(o) - 4M - <Mo) + 4m] [tPM - 4(V 
^ic(o) + 4M - ^c(o) - 4m] K(o) + 4( 
4M4 2c (o) + V> 2c (o>Mo)" • 



(B9a) 

(B9b) 

(B9c) 
(B9d) 

(BIO) 
(Bll) 

(B12) 
(B13) 



The fermion field i/'s is decoupled from the other fields, and thus we neglect it in the following. We bosonize the 
remaining three fermions again: 



V'icO) = 

^2c( X ) = 
^si x ) = 



D 



2ttUvf 

D 
2nhvF 

D 
2-kTivf 



1/2 



1/2 



1/2 



(B14a) 
(B14b) 
(B14c) 



where the tp's are bosonic fields and the 77's are Majorana fermions. Equations ( B10| )-( B13 ) are then rewritten as 
H ° = ^J {[ V ^W1 2 + I V ^)] 2 + [V^(^)] 2 } dx, (B15) 



Hi 



c 



E c 



Va(0)-Vi(0)-7r \N- 



-i 2 



(B16) 



\/2 V2 
Hv = D\r L \ sin ^(0) - r? 2 sin^ 2 (0)] sinv3 s (0) + D\r R \ cos^O) - tj 2 cos <p 2 (0)] cos <f s (0), (B17) 

I = -i—r\ Y r} 2 sin[^i(0) + y 2 (0)], 



(B18) 



where rjj = irj lc rj s and r/ 2 — iVicVs- ^ n ^1- (B16) the parameter iV is shifted by | so that the Hamiltonian (B15)-(B17) 
yields the same perturbation series for the partition function in powers of |r| as the orig inal Hamiltonian flBl|) (pj|). 

Let us consider the case where N is half- integer (on resonance) and \tl\ = \rR\ = ro (symmetric barriers). At low 
temperatures the difference ^2(0) — fi(0) is almost fixed to be tt(N — i) due to the charging energy. In the strong- 
tunneling limit (ro < 1), we can thus take average of the Hamiltonian over ip 2 — (fi to get an effective Hamiltonian 
for ip s and if = ((ft + ip 2 )j\f2: 
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(//„ - Ih ^ I {[V<^)] 2 + [V<p s (x)?}dx 1 2 ' ' E( 



4ir 
eD 

-i-r-r]^ sin 



7r V irD 



1/4 



g(0) 
x/2 



V^^(0) 



(B19) 
(B20) 



Introducing new bosonic fields, (p a = \j\ i f> — \j\ L P s an d ~ ' 
current as 



§</? s , we rewrite the Hamiltonian and the 



(Ho + Hv) = ^J {[V^(a;)] 2 + [V^(*)] 2 } + ^ ( 



2V2 (lE c ^ 1/A 



irD 



Dr^Gx cos 



1 = a, sin 

ft 



^.(0) + V^(0) 



(B21) 
(B22) 



are Pauli matrices. The Hamiltonian ( B21 ) is the same as that of a single 

3 



where a x = [r\ x — rj 2 )/^/2 and a z — —ir\x} 

impurity in the g = | Luttinger liquidO Since this g = | problem is related to the 4-channel Kondo problem in 
the Toulouse limit ,E£I the above transformations show that in the strong-tunneling limit our bosonized Hamiltonian is 
indeed equivalent to the 4-channel Kondo pr oblem . The infrared stable fixed point of the 4-channel Kondo problem 
corresponds to the case with ro — > oo in Eq. ( |B2l| ). At this point a x cos[(3/2) 1 / 2 (p b ] takes —1 so that the field (p b (0) is 
pinned at either <p b (0) = ±2n(2/3) 1 / 2 7r (for a x = —1) or </? b (0) — 7r = ±2 m(2/3) 1 / 2 7r (for a x — 1), where m and n are 



integers. Thus in calculating the dimension of the current operator ( B22 ) we may neglect </? fc (0). We then immediately 
see,that the dimension is d = 1 because the operator exp[z(4/3) 1 / 2 (^ a (0)] has dimension 2/3 and the dimension of a z 
isL3 1/3. This means that the leading term of the conductance is temperature- independent at T — > 0. To calculate 
the conductance we need to know the coefficient A of the correlator ([I (t) , I (0)]) = —Ae 2 /t 2 at T — > 0, which is a 
difficult problem. We expect, however, that A should be universal, implying that at T — the linear conductance G 
has a universal value wXe 2 /2h independent of the initial value of ro- 
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